Polarized polariton condensates and coupled XY models 
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Microcavity polaritons, which at low temperatures can condense to a macroscopic coherent state, 
possess a polarization degree of freedom. This article discusses the phase diagram of such a system, 
showing the boundaries between differently polarized condensates. The Bogoliubov approximation 
is shown to have problems in describing the transition between differently polarized phases; the 
Hartree-Fock-Popov approximation performs better, and compares well to exact results that can be 
used in the limit where the left- and right-circular polarization states decouple. The effect on the 
phase boundary of various symmetry breaking terms present in real microcavities are also considered. 

PACS numbers: 71.36.-l-c, 03.75.Mn 



I. INTRODUCTION 

The recent experimental progress in attaining sponta- 
neous coherence in a thermalized degenerate gas of mi- 
crocavity polaritonsii^i^iii^i^ has extended the range of 
systems in which quantum condensation may be studied. 
As well as sharing many features with previous examples 
of quantum condensates (such as superfluid Helium, cold 
atoms and superconductors) polariton condensates pos- 
sess naturally a number of distinguishing features (see 
e.g. Refs. and Refs. therein). 

This article considers the combined effect of two such 
features of condensed polaritons: their polarization de- 
gree of freedom and confinement to two dimensions. Re- 
cent worksiSiiiii^ have considered some effects arising 
from the polarization degree of freedom; including a ten- 
tative phase diagram^ ^ derived from a zero-temperature 
mean-field theory, this phase diagram implied a transi- 
tion temperature which vanished at a critical magnetic 
field. This article has two aims: first to examine more 
carefully the phase diagram (critical temperature vs mag- 
netic field) of a polarized polariton condensate, this is 
found to be rather different to that in Ref. [ll|; second 
to discuss how the phase diagram would be affected by 
terms breaking polarization rotation symmetry that are 
expected in real microcavities. 

In order to discuss the effects of polarization on the po- 
lariton phase diagram in the most transparent way, this 
article makes various simplifications: it considers an infi- 
nite two-dimensional polariton system, in thermal equi- 
librium. In addition, the model described in Refs. Hllfl^ is 
used, which is applicable in the limit of very low densities 
and temperatures, where only the low energy part of the 
lower polariton dispersion — with a quadratic dispersion 
— is thermally populated. At higher densities, one must 
take account of the non-quadratic dispersion of lower po- 
laritons, and the possibility of other excitations depleting 
the condensate*. These considerations will change the 
dependence of critical temperature on density, but the 
nature of the possible polarized phases and topology of 
the phase diagram as a function of magnetic field should 
not significantly change. 

Current experiments remain some distance from this 



limit of low temperature, low density, infinite, clean equi- 
librium systems: The experimental densities are at the 
point where other excitations start to become relevanl4, 
the polariton clouds are relatively small and their pro- 
files strongly affected by photonic disorder— i^i^ or trap- 
ping potentials^, and pumping and decay have notice- 
able effects^'^'^. It does however remain a useful exercise 
to understand how spin degrees of freedom modify the 
phase diagram in the ideal system first, before including 
these various extra complications, both to understand 
the nature of the possible phases and transitions, and 
also to have a basis from which one can investigate the 
differences introduced by pumping, decay, disorder and 
finite sizes. In addition, improvements in fabrication of 
microcavity samples can be expected to reduce the ef- 
fect of disorder, and increase the lifetime of polaritons, 
opening the possibility that future generations of experi- 
ments might come closer to the idealized model discussed 
here. A review of how some of these effects of disorder, 
pumping, decay etc affect coherence in current micro- 
cavity polariton experiments can be found for example 
in Ref. M. 



II. THE MODEL 

The polarization of the polariton can be written as a 
two-component complex spinor -0; 
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Here I and r are complex coefficients, describing the state 
in the basis of left- and right-circular polarizations. The 
model of Rubo et al.^ in this basis is: 
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The term ft describes a magnetic field that favors ei- 
ther left- or right-circular polarization. From this point 
onwards h — I. In 2D the form of the phase di- 
agram is controlled by two dimensionless parameters; 
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mUo and {Uq — 2Ui)/Uo. For mC/o, estimates includ- 
ing effects of excitonic disorder— give Uq — 3/ieV(/ini)^, 
and the polariton mass 1/rn = 7600/zeV(/im)^, lead to 
mUo ~ 4 X 10~^. Due to the tendency toward biex- 
citon formation, opposite polarizations of polaritons at- 
tract, hence > Uq — 2Ui > —Uq; at = this im- 
plies |Z| = \r\ describing linear polarization of light^^. 
For the typical values relevant for microcavity polaritons 
\Uo — 2Ui\ ^ Uq; i.e. the interaction between left- and 
right-circularly polarized light is relatively smalU^; a typ- 
ical estimate is {Uq — 2Ui) = —O.IUq. The model of 
Eq. (21) has been studied in the context of atomic con- 
densation, e.g. Ref. [13, where topological defects and 
phase separation were investigated; as discussed there, 
phase separation requires Ui < 0, so is not considered 
here. 



III. CALCULATING CRITICAL 
TEMPERATURE 

In an infinite two-dimensional system, for a single com- 
ponent Bose gas, it is well known that the phase transi- 
tion is a Berezhinskii-Kosterlitz-Thouless (BKT) transi- 
tion. Above the BKT transition vortices proliferate, lead- 
ing to exponentially decaying correlations; below there 
is a non-zero superfluid density, and a quasi-condensate 
density, but at long distances phase fluctuations lead to 
power law decay of correlations. To see how this is mod- 
ified in a spinor Bose gas, one needs to consider the el- 
ementary vortices of the model in Eq. As discussed 
by Rubo^, these elementary vortices are separate vor- 
tices of left- and right-circularly polarized light. Thus, 
to a first approximation, the critical temperature occurs 
when the phase stiffness of Z or of r becomes small enough 
that vortices of / or r will proliferate. Note however that 
although the elementary vortices are separate vortices of 
I and r, as long as Uq — 2Ui ^ then the long wavelength 
phase fluctuations are mixture a^°i^^ of I and r. 

When I and r vortices are independent, the critical 
temperature at which vortices proliferate is given by 
p^f = {2/Tr)mkBT where p^''' are the phase stiffness 
of each species. Section IIVI will show that in the ab- 
sence of symmetry breaking terms beyond the model in 
Eq. ([2]), it is appropriate to treat these transitions in- 
dependently, and will discuss how symmetry breaking 
terms modify this argument. This article will present 
the phase diagram of critical temperature as a function 
of magnetic field, working at a fixed total density of po- 
laritons. With such an approach, to calculate the phase 
boundary one will require expressions for the total den- 
sity /3totai(/^, ^jT) of the coupled system, and for the su- 
perfluid densities for each separate polarization p^'*"; the 
following sections will discuss calculating these quanti- 
ties. 



A. Decoupled case, Uo = 2Ui 

The interaction between left and right polarizations is 
relatively weak, and so it is worth first considering the 
special case Uq — 2Ui, for which the polarizations de- 
couple. This case is simple for two reasons: Firstly, it 
is clear that there are two completely independent phase 
transitions associated with BKT transitions for the / and 
r polarizations. Secondly, one may reuse the equation 
of state pi(/i, T) and critical chemical potential fici for a 
one-component 2D Bose gas, which allows a comparison 
between the exact equation of state and various pertur- 
bative approximations. 
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FIG. 1: (Color online) Equation of state, p(/i,r) = Tf{fi,T) 
for a 2D single component Bose gas, as given by HFP ap- 
proximation and Monte Carlo results of Ref. The wider 
spaced dots for the lower branch of the HFP calculation are 
unphysical; they lead to the corresponding wider spaced dots 
in Fig. H 

The equation of state for a single component Bose gas, 
calculated using a Monte Carlo (MC) method, are de- 
scribed in Ref. I la. In terms of the single component equa- 
tion of state, the density for the two component model 
may be written: ptotai(p, = pi(p - ^,T) + pi(/i + 

ri, r), and so the critical temperatures at field fl are the 
solutions of the equation: 

Ptotai = pi(pci,T'c) + Pi(Pci ± 2n,Tc). (3) 

In two dimensions, the equation of state can be written 
pi(p,r) = Tpi(fj,/T) and the critical chemical potential 
scales as pci = XcT. Given pi{x) (as shown in Fig. [T|) 
and Xc, it is straightforward to find Tc(f2) for the two 
component case: The two-component equation, Eq. ([3]), 
becomes 

Ptotai - [pi{xc) + pi{xc ± 2n/T,)] EE T,F{n/T,). (4) 

The phase boundary calculated from this equation of 
state is shown by the solid line in Fig. [51 The solution as 
Tc — > requires pi — > cx), which for the MC calculation 
occurs at a finite O as ri/Tc — > oo. 

Physically, the phases marked I and r correspond to 
pure circular polarizations, while the phase l + r will have 
an elliptical polarization in general, and at = 0, where 
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the densities of I and r match, it will be linear. Ref. 11 
showed that at T = 0, at the critical fl separating I from 
I + r, there is a quadratic gapless mode corresponding 
to fluctuations of r. It was suggested there that this 
mode prevents superfluidity, and so suppresses the tran- 
sition temperature to zero at this point. However, this 
mode corresponds to excitations of r which have decou- 
pled from those of I at this point, and so a force acting 
only on I can still have a superfluid response; this is ob- 
viously the case when Uq — 2[/i, but remains true even 
for non-zero interactions. 
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FIG. 2: (Color online) Critical temperature vs magnetic field 
(both measured in units of fcsTdcg = no/m), at fixed total 
density. Labels /, r, l-\-r and neither state which polarizations 
are superfluid. Red solid line: Monte Carlo results valid when 
IJ\ — 0.5(7o- Blue dotted: HFP results for the same param- 
eters. Green dashed: HFP for a weak inter-polarization in- 
teraction. The HFP lines at intermediate temperatures, with 
wider spaced dots and dashes, are unphysical. 
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General case, Uo 7^ 2Ui 



For Uq ^ 2Ui, one would either need to perform new 
numerical simulations, or find an appropriate perturba- 
tion scheme. The simplest scheme one might consider is 
the Bogoliubov approach, which considers only the self 
energy due to interactions with the condensate. This 
approach is however incapable of describing transitions 
such as that between l + r both condensed and just I con- 
densed. In the phase with just / condensed, as /i il, 
the density of r will diverge, and so the density of I (and 
hence the critical temperature) will go to zero. Consider- 
ing this Bogoliubov approximation as applied to the case 
of Uq = 2Ui, this zero of the critical temperature can be 
understood as follows: If one plots the Bogoliubov ap- 
proximation to the equation of state pi{x), there is an 
unphysical divergence as a; = fi/T — > 0~, while the exact 
equation of state is smooth. It is clear that an unphysi- 
cal divergence of pi{x) at a finite value of x leads to an 
unphysical extra solution of Eq. Q at f7 = 0, T = 0. 
This would produce a phase diagram like that of Ref. [Ill, 



but is an artefact of neglecting the self energy due to the 
non-condensed particles. 

A better approach is the Hartree-Fock-Popov (HFP) 
method, see e.g. Refs. [Mil which includes self ener- 
gies due to the population of thermal and quantum fluc- 
tuations, but neglects anomalous correlations {iH'') etc. 
By including self energies due to fluctuation populations, 
there is no divergence of density, as a large density 
would lead to a large self energy, reducing the effective 
chemical chemical potential for r. 

The HFP method divides each polarization into a 
quasi-condensate density, p^^ and fluctuation density 
p^j^ . Both densities contribute to the gap equation: 

Uoip'o + V/) + {Uo - 2Ui)[pl + P}) = P + ^ (5) 
C/o(pS + 2p}) + {Uo - 2UM + p'f)=p-n (6) 

As in Ref. 17, the fluctuation density is found from 
the correlation function ^^^^(x)Z(O)) = p'(x) exp[— A'(x)]. 
The exponent A'(x) describes the power law decay of 
correlations at long distances, and the prefactor p'(x) 
describes decay from -I- at x = to the quasi- 
condensate density Pq at intermediate distances. Thus 
P/ — — p'(oo)- The correlation function is found 

from the the effective action for density and phase fluc- 
tuations. Writing / = \/ Po + Tr'e*^ (and similarly for r), 
one may define — (6*', tt', 6''', Tr*"), in terms of which: 



,55 = ^ (5«'t(w, k)g-\uj, fc)(5*(w, k) 
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Because of the spinor structure, the current-current re- 
sponse functions (and hence the superfiuid density) have 
a tensor structure. For example, left superfluid density 
is given by p[ = {pq + p^j:) — where = mxx is the 
normal density, found from the transverse part of the 
current-current response function: 



[g{k,u;hnk,k)g{k,u;)j^ik,k)] , (8) 



and 7", the current vertex is (for a — I): 



ll{k,k) 



a'2 




0'2 = 



(9) 



The above method gives the total density and super- 
fluid densities when both polarizations are condensed, 
and so gives the lower temperature boundaries shown in 
Fig. [2l The HFP method can also be adapted to find the 
higher temperature boundary between a single polariza- 
tion condensate and the normal state. Consider the case 
when only I is condensed; in this case the total density 



4 



is Ptotai = Po + Pf + P/' ^^"^ t'^^ S^P equation is just 
Eq. ([5]) with Pq = 0. With a single condensate in the 
HFP approximation, there is no coupHng between the 
fluctuations of the I and r polarizations, but only mean- 
field shifts from the quasi-condensate density, hence, the 
inverse Green's function for fluctuations of I is given by 
the top left 2x2 block of Eq. ([7]). For fluctuations of r, 
one has p^ = nsi^k + S), where S incorporates both 
self energy and chemical potential, and is given by 

E = 2Uop} + (L/o - 2UM + P\) - (m - 

Eliminating p using the gap equation, Eq. (O, this be- 
comes: 

S = (C/o + 2C/i)p^ - 2C/ip[, - (C/o + 2Ui)p\ + 2n. 

Figure [5] clearly shows that the HFP approximation 
cannot accurately describe the entire phase boundary, 
even when Uq = 2Ui. Although the low temperature 
boundary calculated by the HFP and exact equation of 
state match well, the high temperature boundaries match 
less well, and worse yet the critical temperatures do not 
join at n = 0. This failure is because the HFP equation 
of state is discontinuous at the transition, and so is in- 
accurate when the high temperature transition is close 
to the critical region of the low temperature boundary. 
However, far from the critical region the HFP approx- 
imation is reasonable, i.e. at least as good as the HFP 
approximation would be for a single component gas. 

Despite its failings, the HFP method is valuable; flrstly 
it shows that the effect of weak inter-polarization inter- 
actions on the critical temperature is small. Secondly it 
reveals an intriguing property of the transverse current- 
current response function at zero temperature. In a single 
component condensate this vanishes, meaning the entire 
system is superfluid. For coupled polarizations it does 
not vanish, but instead one has Xt ~ Xt^ ~ ~Xt ~ 
~Xt 7^ ^- This identity ensures that a force acting on 
the total density has a longitudinal (i.e. superfluid) re- 
sponse, but a force acting on only one polarization need 
not. This is expected, since there is overall Galilean in- 
variance for a change of velocity of both polarizations, 
but not under changes of the velocity of just one po- 
larization. This zero-temperature transverse response is 
however small, mx"/ptotai oc mllo, as it arises due to the 
quantum condensate depletion. 



IV. SYMMETRY BREAKING EFFECTS IN 
NON-IDEAL CAVITIES 

The discussion so far is based on separate BKT tran- 
sitions associated with the proliferation of each kind of 
vortex; i.e. the effective action is that of an XY model: 

fc^ = - E cos(^l - 0^.) + cos{ei - 0^)] . (10) 



where K^'"^ = p^f / (mkBT). This section discusses some 
possible corrections to this effective action that may 
change the critical behavior. Three effects that are con- 
sidered in detail: The first is short range attraction be- 
tween vortices of 6*' and 9'^ due to density-density interac- 
tions, which are already present in the model of Eq. ([2|). 
The other two effects concern reductions in symmetry 
present in real cavities. Even in an ideal quantum well, 
the symmetry group for Zinc-Blende structures is not 
cylindrical, but D2dr^', this means there are a preferred 
pair of axes, and so interactions between polarizations of 
light do not have complete rotation symmetryii. When 
there is asymmetry between the quantum well interfaces, 
this symmetry is yet further reduced to C2ir^^, leading 
to a preferred linear polarization, causing a splitting of 
the quadratic polarization terms. Such a splitting can 
be induced by applying an electric field along the growth 
direction ^^; even without an applied field, such a reduc- 
tion of symmetry is observed in current experimentsii^. 
These two types of symmetry reduction lead to perturba- 
tions that prefer certain phase relationships between the 
I and r fields, and so modify the order parameter space, 
hence change the behavior of the phase transitions. An- 
other type of symmetry reduction concerns splitting be- 
tween transverse electric (TE) and transverse magnetic 
(TM) modes due both to the cavity and exciton-photon 
coupling*-^^'^'^; this leads to an interesting coupling be- 
tween gradients of I and r fields^!, which may shift the 
transition but does not change the symmetry of the or- 
der parameter space so is not considered in detail here. 
The next section describes how the above considerations 
should be incorporated as perturbations to Eq. (fTO|) , and 
Sec. IIVBI then discusses their effect on the phase bound- 
ary. 



A. Nature of perturbations to Eq. (|10p 

The short-range attraction between opposite vortices 
arises because when Uq ^ 2Ui, a vortex of one polariza- 
tion is associated with a density modulation of the other 
(see Ref . [l^ . Thus a configuration with a vortex in each 
polarization has a lower energy when these vortices are 
colocated, independent of whether the phase windings of 
each polarization are aligned or anti-aligned. This energy 
difference AE is finite, and numerical analysis for = 
and smaU {2Ui - Uo)/Uo gives AE/Ec oc {2Ui - Uo)/Uo, 
where Ec is the core energy of a single vortex. This effect 
can be represented in the Kosterlitz-Thouless scenario by 
ascribing a larger fugacity to hybrid vortices (where both 
I and r wind) than to single vortices. 

The reduction of symmetry in real crystals leads to 
terms that couple the phase of I and r fields. Breaking 
the symmetry to D2d means there are two orthogonal pre- 
ferred directions of polarization; since these orthogonal 
directions are equivalent, this does not produce splitting 
in the quadratic terms, but as discussed in Ref. '11;, there 
is splitting due to interactions that may be written as 
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\ipx\'^ + IV'yl^; using Eq. ([T]) this is proportional to 

|/ + r|4 + \l- r|4 = 2{\l\^ + \rff + 4(rr + lr*)\ 

Breaking the symmetry yet further to means favoring 
a specific hnear polarization, and this will introduce a 
splitting at quadratic order; such a term would take the 
form {l*re^^° + H.c). Both of the above terms can be 
written in the notation in Eq. (|10p as a term 

<5i?p-Ap^cosb(0l~0[)], 

i 

where p = 2 describes the reduction to D2d and p = 1 
the reduction to C2V 



one transition. However, numerical simulations^! of a 
closely related model^^ suggest that at if' = there 
are two close but separate transitions: A higher temper- 
ature Ising transition where 9^ — 6^ becomes locked at 





neither ^ 




A 


\ / ppl mised ' \ / 




B 






c 






r 


/ '' / l+r \ ^\ \ 




Magnetic field, Q. 



B. Effect of perturbations on phase diagram 

The generalization of Eq. (fTO|) , including the perturba- 
tion 5Hp was first discussed by Granato et al.-^^. They in- 
cluded hybrid vortices, but only with aligned phase wind- 
ings. Their discussion is based on studying the renormal- 
ization group (RG) flow in the Coulomb gas formulation 
of the model, where the dynamic variables are the po- 
sitions of vortices. Using their formalism it is straight- 
forward to show that when Ap = 0, the effect of hybrid 
vortices is unimportant: While in principle hybrid vor- 
tices can generate long-range interactions between left 
and right vortices, this does not occur if the energy of 
aligned and anti-aligned hybrid vortices are equal, which 
is the case here^^. This means that if Ap = the transi- 
tion is exactly the scenario assumed above, where a BKT 
transition occurs if either species of single vortex prolif- 
erates. 

As the hybrid vortices alone have no significant effect, 
the model is exactly that discussed by Granato et al.— . 
Unfortunately perturbative RG cannot adequately de- 
scribe the model with 5Hp for p < 4; the RG calcu- 
lation always tends to a phase in which some species 
of vortex or dual vortex proliferate, and perturbation 
theory breaks down2£. In addition, since one possible 
scenario discussed below involves two closely separated 
phase transitions, numerical approaches are challenging, 
as large system sizes are needed to prevent diverging cor- 
relation lengths near one phase boundary masking effects 
of the other—. For this reason, rather than a definite con- 
clusion, the remainder of this article discusses the various 
possible scenarios (Fig. [3]) that have been proposed for 
5H2, and their consequence for the polarized condensate. 
Granato et al.= suggest two possible topologies of the 
phase diagram: either the same as without 5H2 (option 
A), or with a region for ~ where there is a di- 
rect transition from the f + g condensate to the uncon- 
densed state (option B). Their suggestion was based on 
an argument^^ that when K'' ~ there can be only 



FIG. 3: (Color online) Possible topologies of phase boundary 
in the presence of p = 2 (Le. quantum well) symmetry break- 
ing; A and B reflect Ref. [25|, C contains the possible polarized 
but non-superfluid phase suggested by Ref. [271 . 

either or tt, and a lower temperature BKT transition 
where power law correlations of average phase occur. For 
polaritons, this scenario (option C) implies an interme- 
diate phase with linear polarization but no superfluidity. 
While polariton experiments might discriminate between 
these scenarios, this would be made difficult by the close 
spacing between the transitions, and the need for large 
homogeneous systems to avoid finite size effects. 

V. CONCLUSIONS 

The phase diagram of polarized polariton condensate 
contains transitions between regions of circularly polar- 
ized condensates and elliptical polarizations. Comparing 
numerical and perturbative methods when the two polar- 
izations decouple reveal limitations for the perturbative 
methods near the critical point, but they suggest that 
weak interactions between opposite polarizations have a 
small effect on the phase boundary. Terms which break 
the symmetry between different linear polarizations can 
change the topology of the phase diagram; such effects 
may however be very small, but present the possibility of 
using spinor condensates for experimental investigation 
of the topology of this phase boundary. 
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